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Abstract 

The Separation of Variables theory for the Hamilton-Jacobi equa- 
tion is 'by definition' related to the use of special kinds of coordinates, 
for example Jacobi coordinates on the ellipsoid or Stackel systems □ 
in the Euclidean space. However, it is possible and useful to develop 
this theory in a coordinate-independent way: this is the Geometri- 
cal Theory of Separation of Variables. It involves geometrical objects 
(like special submanifolds and foliations) as well as special vector and 
tensor fields like Killing vectors and Killing two-tensors (i.e. isome- 
tries of order one and two), and their conformal extensions; quadratic 
first integrals are associated with the Killing two-tensors. In the re- 
cent years Separable Systems provide mathematical structures studied 
from different points of view. We present here a short review of some 
of these structures and of their applications with particular consider- 
ation to the underlying geometry. Algebraic relations among Killing 
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tensors, quadratic first integrals or their associated second order dif- 
ferential operators and some aspects of approximation with separable 
systems are considered. This paper has been presented as a poster at 
Dynamics Days Europe 2008, Delft 25-29 August 2008. 

Introduction 

There are two levels of the Geometrical Theory of Separation of Variables 
(GTSOV): (i) general or non-orthogonal separation (i.e., not necessarily or- 
thogonal), (ii) orthogonal separation. A geometrical characterization (in 
terms of foliations and Killingvectors and tensors) of the general separation 
has been proposed by Benenti □ . Among the applications of such a character- 
ization we find a coordinate-independent proof of a theorem of Kalnins and 
Miller 0: on a Riemannian manifold with constant curvature the SOV always 
occurs in orthogonal coordinates. A second application is a finer classifica- 
tion of the orthogonal separable systems. Eisenhart's theorem on Stackel 
systems provides a first result on the geometrical characterization of the 
Orthogonal-SOV (OSOV). However, the necessary and sufficient conditions 
written in the Eisenhart statements are redundant. Minimal conditions have 
been proposed by Benenti. According to one of the various possibilities, 
OSOV of the geodesic flow on a Riemannian manifold occurs if and only if 
there exists a Killing tensor K with simple eigenvalues and normal eigenvec- 
tors. Such a tensor has been called a characteristic tensor. Furthermore, if 
a potential energy V is present, it must satisfy the characteristic equation 
dK dV = 0. With respect to Cartesian coordinates on an Euclidean space 
this equation reduces to the so-called Bertrand-Darboux equation. The 
existence of a characteristic KT implies the existence of other n—1 linearly in- 
dependent KT with common eigenvectors which are associated with n — 1 
quadratic first integrals of H by Hi = ^Kj l pjpi + Vi with dVi = K t dV. 
In many cases (for example when the Ricci tensor is null) it is possible to 
build quantum symmetry operators for the Schrodinger equation asso- 
ciated with H by putting Hi = V ' j{KfVi) + Vi and obtaining multiplicative 
separation for the Schrodinger equation in the same coordinates as for the 

3 S. Benenti, J.Math.Phys. 38 6578 (1997) 
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classical system 0. In the following, some few topics of GTSOV are sketched, 
together with some of their possible applications. 

1 Benenti systems 

Among the orthogonal separable systems we find the special class of the 
so-called Benenti systems (or L-systems) 0, for which the separation is 
characterized by a special conformal Killing tensor L called Benenti tensor 
(or L-tensor). For these systems a complete set of quadratic first integrals 
can be constructed by a pure algebraic and coordinate independent method 
starting from L. The geodesic flow on an asymmetric ellipsoid (Jacobi) is 
an example of a Benenti system. Several recent papers have shown that 
Benenti systems have a very rich structure closely related to other fields of 
research. An L-tensor is a conformal KT of order (1,1) which is torsionless 
with (real) simple eigenvalues. Let L a symmetric 2-tensor, then the tensors 
K a , a — 0, . . . , n — 1, defined by 

K = I, K a = -tr(K a _ 1 L)I-K a _ 1 L, a > 1, 
a 

are n independent KT with common normal eigenvectors iif L is a L-tensor. 
The KT's K a determine an orthogonal separable system, a Stackel system 
called Benenti system. Not all the Stackel systems are Benenti systems. The 
first integrals are obtained from the KT as exposed in the Introduction; a 
recursive formula similar to the previous one can produce the potentials in 
each of the first integrals. 

1.1 Separable coordinates for triangular Newton equa- 
tions 

(see 0) Two dynamical systems on the same configuration manifold Q are 
equivalent if their motions on Q locally coincide as unparametrized curves; 
i.e. up to time transformations of the kind dt = fi{q)dt' ( 7 ). For natural 
systems: a dynamical system admits a Lagrangian equivalent system iif the 

6 S.Benenti, C.Chanu and G. Rastelli, J.Math. Phys. 43 5183 (2002) 
7 S. Benenti, Acta Applicanda Mathematica 87 33 (2005) 

8 K.Marciniak and S. Rauch-Wojciechowski, Studies in Applied Mathematics 118 45 
(2007) 
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force F is such that F = — A~ 1 W, where A is the cofactor of a special 
conformal Killing tensor J: A = (detJ)J^ 1 . A special conformal KT with 
simple eigenvalues is a L tensor. An L tensor is a special conformal KT. A 
cofactor pair system is a cofactor system in two distinct ways. A cofactor pair 
system such that J has real simple eigenvalues is equivalent to an L system 
generated by J. The triangular Newton equations 

qi = M i (q 1 ,...,q i ) i = l...n 

are not in general a Lagrangian system, but if M = —{cof G)~ 1 'Vk, where 
G is the matrix associated with the characteristic KT of elliptic coordinates 
in the Euclidean space, it is equivalent to a Lagrangian one. In this case in 
fact G is an L tensor, then a special conformal KT and generates a Benenti 
system. The system equivalent to the TNE is then Hamiltonian and admits 
SOV in the orthogonal coordinates of the Benenti system. 
Example : 

Let (xi,x 2 ) Cartesian coordinates, let the triangular system be X\ = 
— Ax i, X2 = Qx\ — 4x2] this is a cofactor system with G = ( ^ 2x^j ' 

The separable coordinates (not a Stackel system of the Euclidean plane) are 
Ui = Xi and u 2 = \x\ + X2 and their associated quadratic first integrals 
provide the quadrature of the system. 

1.2 Separation curves, Stackel systems and soliton sys- 
tems 

(see 0) A set of n relations of the form (f>i(Xi, fa, a\, ■ ■ ■ , a n ) = dj 6 
M, such that (Xi, fa) are canonical coordinates and det(<90i/<9aj) 7^ are 
called separation relations, when they are all of the same form are called 
separation curve. If the separation curve is of the form 

H^- 1 + H 2 X n - 2 + ... + H n = ^AV + A fc 

where the Hi are the polynomials Hi = ^/i T 'KiG^fi + with G^ = 
L m G (o) ; G (o) = diag (1/Ai; _ ; 1/An); Ai = - Xj) and where L = 

diag (Ai, . . . , A n ) is a conformal KT with respect to the metrics G^ m \ The 
9 M.Blaszak and K. Marciniak, JMP 47 032904 (2006) 
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tensors K { are given by K i+ i = LKi + qj, K x = I, K n+ i = where the 
<7i(A) are coefficients of the characteristic polynomial of L. The potentials 
are obtained by the recursion Vf = V^ 1 - qiVf' 1 . Therefore, the H t form 
a Benenti system generated by the L-tensor L with separable coordinates 
Ai. The qi are new coordinates on the configuration manifold and for the 
Benenti systems of above a hierarchy of integrable dispersionless equations, 
called Killing dispersionless system, whose solutions are called in general 
solitons, is obtained from the equation -^-qj = Ki(q)-^q where ti is the 
evolution parameter (time) of the Hamiltonian flow of Hi and x — t\. 
Example: 

H 1 X + H 2 = |A/i 2 + A 4 be the separation curve, the diaganalized Hamil- 
tonians (Benenti system) are H 1 = 2 (a 1 -a 2 ) — ^2^2 + 2(A 4 — Af), H 2 = 
2 (A 2 - 2 Al ) (^i - iA + 2 (A? - \ 3 2 ). qi = Ai + A 2 , q\ = -4(AiA 2 ) where (A;, fj) are 
canonical Parabolic coordinates. In Cartesian canonical coordinates (qi,p l ) 
they are 

1 11 11 

H i = g(Pi + P2) + 9i + T^i^, #2 = 2 (92P1P2 - <?iP 2 + g?2 + 2^i^)> 

ifi is the Hamiltonian of a separable case of the Henon-Heiles system. The 
integrable dispersionless PDE are 

1 1 

Qi,t = 2^2?2,x, ?2,t = g 925m - QlQ2,x 

that, after differentiation and elimination of q 2 yield the KdV soliton equa- 
tion q 1>t + \q\ tXXX + 3qiq\, x = together with one differential consequence of 
it. 



2 Superintegrability 

Classical Hamiltonian systems in n dimensional manifolds admit at most 
2n — 1 functionally independent (f.i.) first integrals. They are Liouville inte- 
grable if n f.i. first integrals are in involution. Integrable systems with more 
than n f.i. first integrals are called superintegrable. Well known examples 
are the harmonic oscillator and the Kepler-Coulomb system. Because the in- 
tegral curves of the systems stay in the intersection of the level hypersurfaces 
of the first integrals, the occurrence of the maximum number of functionally 
independent first integrals, i.e 2n — 1, implies that the integral curves of the 
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system can be determined by algebraic methods. Hamiltonian systems sepa- 
rable in multiple coordinate systems admits often more than n f.i. quadratic 
first integrals and they are said to be quadratically superintegrable. The 
Kepler-Coulomb system is an example of quadratically maximally superinte- 
grable system. By searching among multiseparable systems many instances 
of superin tegra bility have been found recently and their classification is in 
progresflmHB 

2.1 Superintegrable 3-body systems on the line 



(see Hj) By using Cylindrical coordinates (r, if), z), with rotational axis z, and 
by indicating with (p r ,Ptp,Pz) their conjugate momenta, let us consider the 
natural Hamiltonian with potential 

The potential V is separable, and therefore Liouville integrable, in Cylin- 
drical, Spherical, Parabolic, Ellipsoidal Prolate and Oblate coordinate sys- 
tems. The five quadratic first integrals associated with the multisepara- 
bility take the form H = § (f r + ±p% + p 2 z ) + ™, H x = \ V % + F(if)), 

H2 = \pl H : ' 



\rp z - zp r f + (l + S) Pi] + (l + S) F(iP), H 4 

\ { z Pr + ^iP\ ~ r PrPz) + ^F(if)). Only four of them are f.i. The first three 
integrals allow the separation of the system in cylindrical coordinates. Let x % 
i = 1 ... 3 the positions of three points on a straight line, (x l ) can be inter- 
preted as Cartesian coordinates of a single point in the space. It is possible 
to show that all interactions among the points of the form 

V = T.jf F '(%^) >=l,.,3(mod3), 



10 E. G. Kalnins, G. Williams, W. Miller Jr and G. S. Pogosyan, J. Math. Phys. 40 708 
(1999) 

n S. Gravel and P.Winternitz, J. Math. Phys. 43 5902 (2002) 
12 C.Daskaloyannis and K.Ypsilantis, J. Math. Phys. 47 042904 (2006) 
13 S. Benenti,C. Chanu and G. Rastelli, J. Math. Phys. 41, 4654 (2000) 
14 P.Tcmpesta et al. eds: Superintegrability in Classical and Quantum Systems, CRM 

proceedings & Lecture Notes 37 Amer. Math. Soc. (2004) 

C. Chanu, L. Degiovanni and G. Rastelli: Superintegrable three body systems on the 

line JMP 49, 112901 (2008) 
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are superintegrable with 4 f.i. first integrals, where Fj are arbitrary functions 
of two variables and X { = x % — x l+1 , i = 1,2,3 (mod 3), r cosip = ~j^( xl ~ 
x 2 ), rsmip = -^(x 1 + x 2 — 2x 3 ), z = -^(x 1 + x 2 + x 3 ). Indeed, it can be 
proved that such V are in the form (pQ), then, admit all the first-integrals Hi 
four of which are always functionally independent. 
Examples: 

The Calogero system 

k\ k 2 k 3 ki 

V ° = ( X l- X 2)2 + (a;2 _ ^3)2 + ( X 3_ x l)2 = ^ Xf ' * & ' 

The Wolfes system 

v kl I h I h 

W (x 1 + x 3 - 2x 2 ) 2 (x 2 + x 1 - 2x 3 ) 2 (z 3 + x 2 - 2X 1 ) 2 

= Yl ki+1 ("X 1 ~ ~JTJ ' * = 1 >---> 3 ( mod 3 )- 

And a new one 

V — — — ki ( ^ + — (mod 3) 

~( X i + X i+1 ~[ X i+2 \ X i+2 X i+2/ 



3 Stackel approximation 

Several attempts have been made to use the Stackel systems as suitable 
integrable Hamiltonians to perturbe in order to approximate given physical 
systems. For example and for the dynamics of stars in elliptical galaxies. 
There, for given separable systems and by using the separable coordinates 
themselves, the possible separable potentials (i.e. satisfying the condition 
d(K dV) = 0) are analyzed with the purpose to find the separable dynamics 
closest to the gravitational potentials of elliptical galaxies. The best fitting 
integrable Hamiltonian is then determined by using perturbative methods. 

16 J. Binncy, S. Trcmainc, Galactic Dynamics, Princeton University press (1994) 
17 P. T. De Zeeuw, D. Lynden-Bell, MNRAS 215 713 (1985) 
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3.1 Decomposition of a potential into integrable and 
perturbative terms 

(see 0) Given a natural (non integrable) Hamiltonian H we search for the 
a separable system (Hi), i = l...n representing "the closest", among all 
separable systems, approximation to the dynamics determined by H. Due to 
the structure of the separable systems, the quantities mj = {H, Li} = dLi/dt 
are linear homogeneous polynomials in the momenta. Then, the scalars /i, = 
gjimlm\ provide informations on the time- variation of the functions Lj along 
the integral curves of H . It is remarkable that d\X{ = d Kid V . We can 
conjecture that the minimization of all these quantities, which are coordinate- 
independent, among all the Stackel systems in some neighborhood of the 
configuration manifold determines the separable system which provides the 
best approximation of the given (nonintegrable) Hamiltonian dynamics in 
the same neighborhood. 
Example: 

Let us consider the Quadrupole field in two dimensions: H = \(j? x + 
+ ^ + ^ — where G, D are real constants. By applying the criteria 
of above we determine the separable Hamiltonian system associated to elliptic 
coordinates with foci of cartesian coordinates, if D/C > 0, x 2 = 2D/G, y = 
and with potential W e = szn (Gv — 4^) , where (u, v) are the same elliptic 
coordinates. The dynamics of this system can be compared with the original 
one and with the system separable in Polar coordinates 18 . 



4 Et cetera 

Several advances in GTSOV cannot take place in this short presentation. 
It follows a short list just to mention them and to provide some references. 
Conformal separatiorfj£@, Stackel transformati 



18 G.Rastelli Decomposition of scalar potentials of natural Hamiltonians into integrable 
and perturbative terms: a naive approach J. Phys.: Conf. Ser. 128 (2008) 012029 (12pp) 
19 S. Benenti, C. Chanu and G. Rastelli, J. Math. Phys. 46 042901 (2005) 
20 M. Chanachowicz, C. Chanu and R. G. McLenaghan, J. Math. Phys. 49 013511 
(2008) 

21 E. G. Kalnins, J. M. Kress and W. Miller Jr, J. Math. Phys 46 053509 (2005) 
22 A. Sergeev and M. Blaszak, J. Phys. A: Math. Theor. 41 105205 (2008) 
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Classification 
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of separable systems in constant curvature manifolds (3-dim Minkowskp, 
2-dim De Sitter and anti De SitteiF*!). 



23 J. T. Horwood and R. G. McLenaghan, J. Math. Phys. 49 023501 (2008) 
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